Tight-binding formulation of the dielectric response in semiconductor nanocrystals 
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We report on a theoretical derivation of the electronic dielectric response of semiconductor 
nanocrystals using a tight-binding framework. Extending to the nanoscale the Hanke and Sham 
approach [Phys. Rev. B 12, 4501 (1975)] developed for bulk semiconductors, we show how local 
field effects can be included in the study of confined systems. A great advantage of this scheme 
is that of being formulated in terms of localized orbitals and thus it requires very few computa- 
tional resources and times. Applications to the optical and screening properties of semiconductor 
nanocrystals are presented here and discussed. Results concerning the absorption cross section, the 
static polarizability and the screening function of InAs (direct gap) and Si (indirect gap) nanocrys- 
tals compare well to both first principles results and experimental data. We also show that the 
present scheme allows us to easily go beyond the continuum dielectric model, based on the Clausius- 
Mossotti equation, which is frequently used to include the nanocrystal surface polarization. Our 
calculations indicate that the continuum dielectric model, used in conjunction with a size depen- 
dent dielectric constant, underestimates the nanocrystal polarizability, leading to exceedingly strong 
surface polarization fields. 
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I. INTRODUCTION 

The role of the external environment on the electronic 
properties of semiconductor nanocrystals has been rais- 
ing a wide interest. This is due to an enhanced con- 
trol of the synthesis techniques that allows for an op- 
timization of the nanocrystal optical performances. In 
particular, it has been shown that silicon nanocrystals 
synthesized in organic solutions have a strong photolu- 
minescence in the blue spectral region, opening the route 
toward biomedical applications of water soluble silicon 
nanoparticleSfii^i^ From the theoretical side, few studies 
have been done on the influence of the external environ- 
ment on the semiconductor nanocrystal properties. In- 
deed, while isolated nanocrystals have been studied in the 
past using the independent particle approximation,i^ii 
the role of local field effects (LFEs) has been analyzed 
within quantum mechanical schemes only recently and 
for small structures^i^ii^ 

A mixture composed of nanocrystals embedded into a 
dielectric background is usually described using contin- 
uum dielectric models (CMs). Such models are based on 
the assumption that a nanocrystal can be approximated 
as a continuum medium with a defined position indepen- 
dent static dielectric constant (it can be either the bulk 
value or a size dependent function), that abruptly goes 
to the background dielectric constant at the nanocrystal 
surface. A basic assumption of CMs is that going from 
the nanocrystal center to the external background, the di- 
electric constant has a sharp discontinuity across the sur- 
face. In the case of non interacting nanoparticles, these 
models reduce to the Maxwell-Garnett equation a^^'^^i^^ 
coupled to the Clausius-Mossotti equation for includ- 
ing the surface polarization fieldsji^ Recent studies have 
shown that CMs are not always adequate to give a fair 
description of small nanoparticles. In fact, by looking 
at the local permittivity profiles, it emerges that close 
to the surface there is a "dielectric dead layer," and 



the matching between the dielectric constant inside the 
structure with that of the background takes place within 
a finite lengtbi^ii^iiiii^ whose value has a primary im- 
portance in the technological applications of semicon- 
ductor nanocrystals. The reduction of the local dielec- 
tric constant near the surface is thus the principal cause 
of the reduction of screening in nanocrystals, whereas 
the contribution due to the opening of the band gap 
has a minor roleii^iiiii^ Moreover, in small and not well 
shaped nanocrystals, both quantum effects^ and atom- 
istic features^S become significant, implying that calcu- 
lations of LFEs within atomistic, quantum mechanical 
frameworks are definitely required. Unfortunately, the 
computational cost of this operation rapidly explodes 
with the nanocrystal size because of the need of com- 
puting a large number of excited states. This is the main 
reason why first principles methods can be used only for 
very small nanocrystals in a size range that is often signif- 
icantly below that of the experiments i^i^ A good com- 
promise between accuracy and the possibility of studying 
nanocrystals comprising thousands of atoms is obtained 
by using a real space approach. In the case of dynamic 
response calculations, time dependent density functional 
schemes based on a real space, real time approach are 
the most promising tools for calculating ab initio the op- 
tical properties of confined systems with a favorable time 
scaling^^^e 

Our choice in this paper is to use a tight-binding 
method for determining the dielectric response of con- 
fined systems. It is important to stress here that the 
tight-binding method must be intended in the widest 
sense and that is the use of a localized basis set. The 
general framework described in Sec. [TT] and the explicit 
expressions derived for both the independent particle po- 
larizability and the real space dielectric function can be 
used with any localized basis set, including the maximally 
localized Wannier functionsi^ Although the explicit ex- 
pressions we derived in this work are very general, the 
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applications we have done to illustrate the theory are lim- 
ited to the simplest version of the tight-binding method 
where the Hamiltonian matrix elements are fitting pa- 
rameters determined by the material bulk band struc- 
ture. We present numerical results for the optical and 
the screening functions of both silicon and indium ar- 
senide nanocrystals. The good agreement between our 
results and the experimental data on one side and the 
results of first principle calculations on the other gives 
an indication that the chosen route for describing the 
dielectric response at the nanoscale is very useful. An 
interesting comparison of the nanocrystal polarizability 
calculated using the full tight-binding approach and the 
one obtained using the Clausius-Mossotti equation gives 
important indications about the range of validity of the 
dielectric medium theories. 



II. GENERAL FRAMEWORK 

In the following, the linear response theory is reviewed 
using a localized orbital basis set. Our derivation is 
based on the framework successfully applied by Hanke 
and Sham to bulk semiconductors.— The main differ- 
ence consists in the fact that the theory is here developed 
for confined systems that, as such, do not have period- 
icity. As a consequence, a different procedure for the 
calculation of the dielectric constant must be employed. 
We will show that, neglecting local fields, the standard 
expression for the nanocrystal dielectric constant within 
the independent particle random phase approximation 
(RPA) is retrieved. ^'^ Instead, an expression is obtained 
when local fields are taken into account. We refer to the 
first approximation (neglect of local fields) as RPA, and 
to the second approximation (inclusion of local fields) as 
RPA-I-LF. Since the main point we wish to address in 
this paper concerns the role of LFEs in both the opti- 
cal response and in the screening, we neglect excitonic 
effects. Besides, while several interesting studies of the 
excitonic effects in nanocrystals have been done in the 
past^Si^i^ systematic studies on the influence of LFEs 
are still lacking. The interest in this field is also moti- 
vated by studies of semiconductor surfaces done in the 
past, which showed the important role played by LFEs 
in that situationi^i^^iM 



A. Tight-binding approach 

As mentioned before, the tight-binding method is a 
powerful tool for the study of confined systems. Al- 
though modern computational facilities and parallel pro- 
gramming lead the calculations to be much more efficient 
than in the past, first principles studies of nanocrystals 
containing several thousand atoms are nowadays an al- 
most impossible task. The situation gets even worst if 
one needs to study the optical response of a nanostruc- 
ture with the inclusion of local fields, since a calculation 



of a large number of excited states has to be performed]^ 
The tight-binding method has all the advantages of a 
formulation in terms of localized atomic orbitals whose 
most important feature is the use of a relatively small ba- 
sis set, which implies a massive sparsity of the relevant 
matrices. With appropriate algorithms for the storage, 
diagonalization, and inversion of a matrix, it is possible 
to implement an extremely efficient computational tool. 



B. Linear response theory 

The starting point is, of course, the diagonalization of 
the single particle Hamiltonian. The nanocrystal wave 
functions are written as a linear combination of localized 
atomic orbitals. 



i^„(r) = ^B,„(R)0,R(r). 



(1) 



crR 



Here, a labels the atomic orbital symmetry and R's are 
the atomic coordinates in the nanostructure. The imple- 
mentation that we use in our calculations is that built 
on the third-nearest-neighbor s'p^ parametrization (with 
(T = s,px,Py,Pz) which has been shown to give a good 
estimation of energy gaps and effective masses of both sil- 
icon and indium arsenide, the two materials that we have 
chosen for illustrating the theory.— i^^i^i^l We neglect the 
spin-orbit interaction and assume that the basis set is 
composed by real functions. In any case, a generalization 
to include the spin-orbit interaction is straightforward. 

The real space independent particle polarizability is 
defined as^ 



P(r,r',c^) 



E 



Ea — huj — ITj 



A^{v)A^{y'). (2) 



The index a ~ {n, n') runs over all the possible tran- 
sitions between the single particle eigenstates. We de- 
fine A„^„,(r) = V«(r)V-'n'(r), £^n,n' = - En, and 
fn.n' = fn' — fn- fn IS the n*'' Icvcl occupation number 
and 7y is a small energy giving rise to a finite broadening 
of the absorption spectra. Using Eq. ([1]), it is a simple 
matter to see that 

P(r, r', u;) = ^ P^,, {^) Ax{r)A,{r'), (3) 

where A — {a, cr', R) labels a pair of atomic orbitals cen- 
tered on the same site.^-i^ In Eq. ([3]), we have intro- 
duced the tight-binding representation of the indepen- 
dent particle polarizability. 



Ea — fioJ — irj 



(4) 



In Eqs. © and (g]), we have defined ^a(i") = 
0CTR(i')0<7'R(r) and Cx,a = S<7n(R)S^'n'(R), respec- 
tively. In order to simplify the calculations, we make 
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a number of approximations that are strictly related to 
the parametrization we have used. The first approxima- 
tion is the neglect of the overlap between atomic orbitals 
localized on different sites. The very small impact of 
this assumption has already been checked in the calcu- 
lation of the optical absorption spectrum of both bulk'*" 
and nanocrystalline silicon.'' In both the cases, the agree- 
ment between the theory and the experimental data was 
very good. The second approximation consists in ne- 
glecting the off-diagonal intrasite contributions due to 
atomic orbitals with a different symmetry (terms with 
a ^ a' m. the A's). On this point, there has been a wide 
discussion in the literature concerning the necessity of 
introducing additional parameters to take into account 
these contributions . "''^i^^^i'^'^i'^^'^^i'^^ Our experience in this 
matter is that, apart from the theoretical problem of 
the breaking of gauge invariance,— using a parametriza- 
tion with many neighbors (i.e., third-nearest-neighbor 
parametrization) makes the off-diagonal intrasite contri- 
butions negligibleii 

The third, most important approximation is in the use 
of a unique function A{y) for each atomic site, calculated 
by averaging all the functions A\{y) at a given site. With 
this approximation, we are neglecting the dependence on 
the atomic symmetry so that the tight-binding labels re- 
duce to A = R. The averaging is simply given by 

1 4 

(T=l 

Our calculations show that these sets of approximations 
do not introduce significant errors in the RPA dielec- 
tric functions. Within RPA-I-LF instead, it can induce 
an error in the calculation of local fields that, accord- 
ing to the quantity one wish to study, can be either ne- 
glected or not. For instance, we shall see below that these 
approximations are responsible of a small difference in 
the screening when our results are compared with den- 
sity functional theory results. Nevertheless, it must be 
stressed that Eq. (O yields a huge simplification, since it 
leads to a significant reduction of the size of the matrix 
to be inverted. In the case of the sp^ parametrization, 
the tight-binding polarizability matrix Px./^ is reduced by 
a factor of 16. 

It is important to underline the role of the Si-H in- 
teraction parameters. All the electronic properties can 
be very sensitive to them, especially in the case of small 
nanocrystals. For instance, in SissHse the gap energy 
increases by more than 10%, with a significant decrease 
of the static dielectric constant and recombination rates, 
when the Si-H interactions change from the values of Ref. 
[35I (obtained by fitting the Sill4 experimental gap) and 
the parameters of Ref. 113 (chosen to give a stronger pas- 
sivation). Throughout this paper, we use the parameters 
of Refs. m and [H, which provide a good agreement of 
the gap energy with other theoretical results and exper- 
imental data. 



C. Dielectric function 

The real space dielectric function in the random phase 
approximation is given by^ 

e(r, r', uj) ^ 5 {r - r') ~ J driu (r, ri) P (ri, r', uj) , (6) 

where u is the bare Coulomb interaction M(r,r') = 
e^/|r' — r| and the integral is done over the whole space. 
Using Eq. ^ for the polarizability, the dielectric func- 
tion reads 

e(r, r', c.) = ^ (r - r') - ^ Jx{r)Px^, (u;) A, (r') , (7) 

A, /J. 

where we have defined the Coulomb integral 

Jxir) = y"driu(r,ri)AA(ri). (8) 

In order to calculate both the optical response and the 
screening function, the inverse of the real space dielectric 
function is required. This kind of calculations can be a 
formidable task if one adopts a first principles point of 
view,^^ Instead, using the present formulation, the inver- 
sion is done with a modest effort. Using simple matrix 
properties, ^^'"^^ the inverse dielectric function can be 
explicitly written as 

e-\r, r', ^) = <5 (r - r') + 5] Ja (r)^^,;. (c.) A^ (r') . (9) 

In this expression, S is the screened polarizability (also 
defined as screening matrix^*) that in the tight-binding 
representation is 

where the dielectric function is given by^ 

i-x^^) ^ ^x^ii-Y^Ux^vPy^uj), (11) 

and the Coulomb interaction matrix is 
Ux,, = J dr Ja (r) A, (r) = J drdv'^^^^^. 

The main point here is that the inversion of the dielec- 
tric function has been reduced to the inversion of the 
dielectric matrix defined in Eq. (|lip with enormous ad- 
vantages with respect to the direct inversion of Eq. ([6]). 
It is remarkable that the derivation of Eq. ([9]) from Eq. 
([S]) is exact, allowing us to reduce the inversion of a large 
matrix to that of a much smaller onci^ The physical mo- 
tivations beyond such a significant reduction of degrees 
of freedom are rooted on the fact that semiconductor 
systems can always be studied using few basis functions 
per atom, as the tight-binding based experience and the 
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Wannier function approaches have widely shown over the 
years. 

Consistent with the approximations outhncd in sec. 
IIIB| we parametrize the matrix elements in Eq. (|12|) 



scopic electric field. In formulas, it is defined ai 



M. 



R,R' 



eV|R'-R| if R^R' 
[/„, if R = R'. 



(13) 



where Uat is an on-site Coulomb interaction term that 
only depends on the atom located at R. In the applica- 
tions of the method discussed below, we use the orbital- 
averaged function AR,(r) defined in Eq. ^ for the cal- 
culation of the on-site Coulomb terms Uat- There is also 
another procedure for the calculation of the on-site pa- 
rameters, often used in literature. It consists in the cal- 
culation of Uat as an average over all the Coulomb inter- 
action energies between pairs of atomic orbitals. ->^i^i^ 
Below, we shall compare the results obtained with the on- 
site terms calculated using either the first or the second 
method. 



III. DIELECTRIC CONSTANT 

The dielectric constant is a macroscopic quantity, 
well defined for extended and periodic systems. How- 
ever, its meaning in the case of a nanocrystal has been 
questioned;^ and new schemes have been proposed for 
the calculation of a local and position dependent dielec- 
tric constant li^iii One of the most interesting results in 
these studies is that well inside a nanocrystal, the local 
dielectric constant is just the same as that of the corre- 
sponding bulk system. However, upon approaching the 
surface, it rapidly changes until matching the value of 
the material supporting the nanocrystal. The immedi- 
ate consequence of this behavior is that the reduction 
of the dielectric constant in a nanocrystal with respect 
to the bulk, observed not only theoretically^'^-iZ, but even 
experimentally,— has to be ascribed to a surface eSec^ 
which has a dominant weight with respect to the band 
gap blueshift contribution.^^ It is therefore evident that 
the surface polarization plays a key role in the dielectric 
response of confined systems. Within a real space view 
of the problem, the surface polarization gives a strong 
discontinuity of the real space polarizability across the 
nanocrystal surface. A key point in this analysis is the 
validation of the continuum models. These models are 
based on the idea of studying a nanocrystal as a uniform 
dielectric sphere. This is a central point when analyzing 
experimental data simply because it is often not clear 
what is the real range of validity of these models, partic- 
ularly when applied to small structures. With our pro- 
cedure, we are able to calculate local fields within a fully 
atomistic quantum mechanical scheme, checking, in this 
way, the CMs upon changing the nanocrystal size. 

The macroscopic dielectric constant eM(w) is defined 
as the response of a system to a long wavelength macro- 



eA/(w) = lim 



1 



q^oe i(q,q,w)' 



(14) 



where the Fourier transform of the inverse dielectric func- 
tion is defined as 



-i(q,q',c.) = l^drdr'e-i 



(r, r',w)e" 



iq r +zqr 



(15) 

We calculate the Fourier transform of real space quan- 
tities, localized on a finite structure, approaching a con- 
stant value outside the structure. In particular, the di- 
electric function defined in Eq. (fT4)) is equal to the dielec- 
tric constant of the background embedding medium in 
the region of space outside the nanocrystal. Here, we con- 
sider spherical nanocrystals in vacuum, so e(r, r',aj) = 1 
when |r|, |r'| > R (i? is the radius, and the coordinate 
system is placed at the center of the nanocrystal). 

In order to make spatial integrations, we introduce an 
integration volume Q. By using values of much greater 
than the nanocrystal volume (taking the limit f2 oo), 
we retrieve Jl-independent results. A main point is that 
the dielectric constant calculated by Eq. (fH)) is not 
just the nanocrystal dielectric constant, since it includes 
the vacuum space encompassing the nanocrystal. Defin- 
ing the filling factor i—ils/i^ as the ratio between the 
nanocrystal volume fls and the integration volume fl, we 
can write the dielectric constant defined in Eq. as 
the average between the nanocrystal dielectric constant 
es and that of the vacuum space {eout — 1)7 

CM = (1 - /) + fes. (16) 

From Eq. (|16p , the nanocrystal dielectric constant 65 can 
be calculated. 

By substitution of Eq. ([6]) into Eq. (jisp. the following 
expression is obtained: 



e ^ (q,q,w) = 1- 



■u(q) 

n 



J2M<DSx,^iio)A^iq), (17) 



where the Coulomb interaction takes the usual form 
w(q) = ATTe^/q^ and A{q) is the Fourier transform of 
A{r). In the q limit, it can be expanded using the 
dipole approximation 



Ax{q) ~ 1 - iq-DA, 



(18) 



where the dipole matrix elements in the tight-binding 
basis set are given by = / dr [a;Q-A;y(r)], a being the 
Cartesian component x, y or z. The dielectric constant 
defined by Eq. ([M]) is thus 



A, /J, 



(19) 



In order to simplify the notation, in the following we 
do not write the explicit dependence on the Cartesian 
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component a. Since in the examples given below we 
only consider spherical nanocrystal of cubic semiconduc- 
tors, the dielectric tensor reduces to a constant, and 
the dependence on a disappears. However, the result 
can be easily generalized to anisotropic systems. In 
the limit r2 — > oo, the previous expression becomes 
tuiyi) = 1 - 47re2^^^i:)AS'A,^(w)£'^/ri, and accord- 
ing to Eq. (fTB)) . the nanocrystal dielectric constant can 
be finally derived as 



^5]/?A5A,^Mi?M. (20) 



As expected, the integration volume simplifies in all ex- 
pressions and only the nanocrystal volume appears in 
the final expression. We point out that this is the dielec- 
tric constant of a confined structure, defined in linear re- 
sponse theory within RPA with LFEs included (RPA-I-LF 
scheme) . If we neglect local field effects (RPA) , we sim- 
ply take the q — > limit of the reciprocal space dielectric 
function, 



RPA 



(to) 



lini e (q, q, uj) 

q-*0 



(21) 



retrieving a standard expression for the RPA dielectric 
constant. 
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FIG. f : (Color online) Absorption spectra of colloidal InAs 
nanocrystals. Present tight-binding calculations (red lines) 
are compared to experimental data of Ref. [s^ (black sym- 
bols). A Gaussian broadening has been used. 



es (c^) = I - ^ ^ i?APA,M (^) . (22) 



The difference between the dielectric function with and 
without LFEs is in the polarizability term that is cither 
screened as in Eq. ([20)) or unscreened as in Eq. (|22|) . We 
refer the reader to Ref. [H for a detailed discussion on 
the physical interpretation of the difference between the 
dielectric constant calculated within both the schemes 
and on the role played by local fields. 



Results 

In this section, the optical properties of Si and InAs 
nanocrystals are illustrated. We recently applied the 
method to the calculation of the Si nanocrystal absorp- 
tion spectra, and good agreement with experiments was 
recoveredi^^ We report in Fig. [T]the absorption spectra of 
InAs nanocrystals calculated with the present RPA-f LF 
tight-binding scheme, and compare them with experi- 
mental data taken from recent measurements on colloidal 
nanocrystals. We use a second-nearest-neighbor tight- 
binding parametrization;^ calculated by fitting the band 
structure and the effective masses of bulk InAs. As Fig. 
[T] shows, the agreement between theory and experiments 
is good. 

In Fig. [2l the static polarizability per volume (left 
panel) and the static polarizability (right panel) of spher- 
ical Si nanocrystals are shown. Classically, the polariz- 
ability of an atomic cluster can be calculated starting 



from the bond polarizabilities and including the dipole- 
dipole interactions in order to take into account local 
field effects.^" Within a RPA4-LF scheme, however, local 
fields are already included in the calculation of the di- 
electric constant, and thus the nanocrystal polarizability 
is directly related to it through the expression 



a=-{es-l) 



(RPA + LF) 



(23) 



The RPA-I-LF results are indicated in Fig. [2] with blue 
cross symbols. An alternative way of calculating the 
polarizability is by modeling the nanocrystal as a con- 
tinuum dielectric system. Within this model, the cal- 
culation of the surface polarization effects leads to the 
Clausius-Mossotti expression.^" This approach is often 
used by modeling the nanocrystal either with the bulk Si 
or a size dependent dielectric constant. In order to check 
whether the continuum model works well for nanocrys- 
tals, we use the RPA static dielectric constant, and cal- 
culate the polarizability through the Clausius-Mossotti 
expression'^ ^ 



3r>s 

47r 



es-1 



(RPA). 



(24) 



The results obtained with this approximation are shown 
by red square symbols. In the limit of large nanocrys- 
tals, the RPA dielectric constant converges to the bulk 
value^^i and we retrieve a large size limit for the polariz- 
ability per volume, indicated with a black dashed line in 
the left panel of Fig. [21 This value has been obtained 
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FIG. 2: (Color online) Spherical nanocrystalline Si polariz- 
ability per volume (left panel) and polarizability (right panel) 
as a function of the nanocrystal size. Red squares symbols, 
RPA with the continuum dielectric model; blue cross sym- 
bols, RPA+LF; dashed line, bulk limit. The solid lines are 
curves obtained by interpolation (see text), with the following 
fitting parameters: RPA, a = 3.5 A and I = 1.18, RPA+LF, 
a = 3.25 A and I = 1.77. The green triangle symbols are the 
TDLDA results obtained in Ref. [U. 



using the experimental bulk Si dielectric constant in Eq. 
([M)) .'^^ With the use of an interpolation formula for the 
size dependent dielectric constant'^. 



e(i?) = 1 + 



1 



1 + (a/RY 



(25) 



a fit of the present tight-binding results is obtained. The 
results are shown as full lines in Fig. [2] It can be 
noted that the RPA results combined with Eq. ([M]) 
underestimate the polarizability per volume, in agree- 
ment with our previous dielectric constant calculationsji^ 
On the right panel of Fig. [2l the tight-binding results 
for the polarizability are compared to a first principles 
time dependent local density approximation (TDLDA) 
calculatiouj^ The agreement for small nanocrystals is 
very good. However, upon increasing the nanocrystal 
size, there is a discrepancy between TDLDA and tight- 
binding results. This is probably due to the overesti- 
mation of the bulk Si dielectric constant within LDA 
which leads to an overestimation of the static polariz- 
ability, with an error that increases with the nanocrystal 
size<«i 

It is worth trying to better understand the left panel of 
Fig. [21 From the figure, it emerges that the nanocrystal 
polarizability per volume is always smaller with respect 
to the large size limit. As the recent literature has shown, 
this decrease of the polarizability for nanosized objects 
is mostly due to an overall decrease of the local polariz- 
ability across the surface, the quantum confinement ef- 
fect giving a relatively small contributiou i^^'^^i^^ It is ex- 
tremely interesting even from an experimental point of 
view trying to understand the origin of the differences 



shown in Fig. [2] between the full calculation and the con- 
tinuum dielectric model. There are at least three sources 
of errors. First, the Clausius-Mossotti equation assumes 
that the local field inside the structure is position in- 
dependent, whereas even classical atomistic calculations 
have shown that the surface bonds see a local field very 
different with respect to the local field seen by the bonds 
placed inside the nanocrystal(2S A second difference is 
that the surface bonds have a different polarizability with 
respect to the bulk Sii^iiii^ Third, the full calculation 
takes into account the quantum confinement effect. 

All of these effects can give opposite contributions to 
the polarizability, and it is not obvious whether the CMs 
either overestimate or underestimate the correct results. 
In light of these considerations, our results appear signif- 
icant since they state that the full results lie somehow in 
the middle between the CMs used with either a size de- 
pendent or the bulk limit static dielectric constant. The 
relative error made in using the CMs is of the order of 
2%-5% in the size range considered here. 

A last significant finding is that the polarizability per 
volume calculated within RPA+LF is almost indepen- 
dent of the choice of the on-site parameter Vat , showing 
a robustness of the method versus Vat in the polarizabil- 
ity calculations. 



IV. POINT CHARGE SCREENING 

The localized orbital formalism described so far is par- 
ticularly suitable for the study of point charge screening. 
The large interest in this problem is clearly related to the 
study of the excitonic interaction energy, in which a cal- 
culation of the full dielectric function is requiredi^ Such 
a calculation is very demanding from a computational 
point of view,^^ and both models and semiempirical tools 
are often used in order to simplify the problemi^ From 
the definition of dielectric function, the screened electron- 
electron potential energy is given by 



w{r,r',u;) = J (i^i<^ ^ (r, ri, tj) w;, (ri, r') . 



(26) 



where the bare Coulomb interaction is ^^(r, r') = e'^ /\r' — 
r|. According to Ref. [2l|, a common way to reduce the 
computational cost of the excitonic energy is by intro- 
ducing a screening function, defined as 



w (r, r', uj) 
Vb (r, r') 



(27) 



Within the present scheme, using the equations previ- 
ously introduced for the dielectric matrix, the following 
expansion is obtained: 



e-i (r, r', c.) = 1 + 1 |r' - r| ^ J;, (r) Sx,^ (c.) (r') . 

(28) 
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The screening function represents the dielectric response 
to a point charge placed at position r. A common ap- 
proximation consists of using a spherical average of the 
screening function calculated at r = even for the re- 
sponse to an off-center point charge. The error made in 
the calculation of self-energies and excitonic energies can 
be often neglect edi^i^i So an effective screening function 
can be defined as e(r, w) = l/e~^(0, r, w), whose spheri- 
cal average is 

1 + ^rJ2x.f, J\ (0) Sx.t, {uj) Jf, (r) 

By using the approximations described in Sec. Ill Bl the 
Coulomb integral is reduced to the spherical average of 
the function J\{r) ~ ,/(r — R). We calculate J using the 
Herman Skillman numerical tables for the atomic wave 
functions. ^'"^ With respect to previous tight-binding cal- 
culations, the present approach allows for a description 
of the screening in the whole real space. The effective 
screening function for SiasHsg is shown in Fig. [S] We 
report the results obtained using different values of the 
on-site Coulomb interaction parameter Uat for Si atoms. 
In one case, we used the value Uat = 8.8 eV (green dot- 
dashed line) , this value being calculated according to the 
prescription given above, using Eq. ([5]). In the second 
case, the results are obtained with Uat = 10.6 eV (blue 
solid line), a value calculated as average between all the 
Coulomb energies between pairs of atomic orbitalsi ^'^^i"^^ 
The present results are in reasonable agreement with 
both ab initio calculations^^ (red dashed line) and self- 
consistent empirical tight-binding method^^ (symbols). 
It is important to remark that the peak position calcu- 
lated here has the same value as in the ab initio calcula- 
tions. However, the present approach slightly underesti- 
mates the screening function in the space region close to 
the surface, and this is likely to be due to the approxi- 
mations made and discussed above. It is worth to point 
out that the behavior of the screening function near the 
surface is independent of Uat, since it only depends on 
the dielectric constant of the material. 

In Fig. m the effective screening function for Sii9iHi48 
is reported. With respect to the smaller SiasHag, there is 
an increase of the screening, that is higher on increasing 
the size, in agreement with recent ab initio calculations.— 
The static screening (blue solid line) is here compared 
with the optical screening, calculated at a frequency 
corresponding to the gap energy, that for Sii9iHi48 is 
Eo = 2.53 eV (red dashed line). The difference with re- 
spect to the static case is in the increase of the maximum, 
although the peak position remains unchanged. This is 
consistent with the fact that the real part of the dielectric 
constant at the gap energy has a larger value than in the 
static limit. 

In Fig. O the effective screening function for a 2.0 
nm size InAs nanocrystal is shown. As said above, this is 
only the electronic contribution to the screening, since we 
consider the atomic nuclei frozen at their bulk positions. 
The on-site Coulomb interaction parameters calculated 




I > \ I \ > \ I \ > \ > L 

2 4 6 8 10 12 
r (a.u.) 



FIG. 3: (Color online) Spherical average of the screening di- 
electric function. Blue solid line, present calculation using an 
on-site parameter for Si Uat = 10.6 eV; green dot-dashed line, 
present calculation using for Si Uat = 8.8 eV; red dashed line, 
ab initio calculation of Ref. [2ll : black circles, self-consistent 
tight-binding calculation of Ref. [TgI . 




5 10 15 20 



r (a.u.) 

FIG. 4: (Color online) Spherical average of the effective 
screening function of Sii9iHi48. Blue solid line, static screen- 
ing {uj — 0); red dashed line, optical screening (calculated at 
the gap energy). 



according to the prescriptions given in Sec. Ill Bl are Uat 
— 9.68 and 6.25 eV for As and In atoms, respectively. At 
variance with Si, in the case of InAs nanocrystals there 
are two possible configurations for each size, depending 
on whether an indium (blue solid line) or an arsenic atom 
(red dashed line) is located at the center of the nanocrys- 
tal. As shown in Fig. [51 both the intensity and the peak 
position assume different values for the two configura- 
tions. However, far from the impurity the screening func- 
tion is exactly the same. This apparently strange result 
can easily be explained using a Thomas-Fermi modeh ^^i^^ 
According to it, the point charge screening is composed 
of two contributions. A classical contribution, close to 
the surface, is due to the surface polarization charge in- 
duced by an external test charge located at the nanocrys- 
tal center. It is well described by classical electrostatics 
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FIG. 5: (Color online) Screening function of a 2.0 nm diameter 
InAs nanocrystal. We studied the configuration with an As 
(In) atom placed in the nanocrystal center, indicated with the 
red dashed (blue solid) line. A zoom of the graph is shown in 
the inset. 

and only depends on the nanocrystal dielectric constant 
and radius. Indeed, it is significant that both the con- 
figurations show the same behavior close to the surface, 
due to the fact that they have the same macroscopic di- 
electric constant. Approaching to the impurity site, the 
deviation from the classical model becomes important 
because of short range interactions. At variance with 
the classical contribution, the behavior in the neighbor- 
hood of the impurity (in particular, the peak position 
and intensity) strongly depends on the local environment 
around the impurity site and this is the cause of the dif- 
ference between the two curves in Fig. [5l Therefore, 
InAs nanocrystals show important bulklike local field ef- 
fects due to the difference of the atomic cores. An in- 
teresting point confirming this view is that the surface 
polarization contribution does not depend on the on-site 
Coulomb term Uat- This parameter influences the ampli- 
tude of the screening function around its maximum, but 
it does not change the classical part close to the surface. 



In the inset the detail of the merging of the two curves is 
shown. The point at which this occurs can be indicatively 
defined as the screening radius of the nanocrystal.— 



V. CONCLUSIONS 

In this paper, a tight-binding framework for the study 
of the dielectric response in confined systems has been 
described. The interest of the method consists in the 
easy inclusion of local field effects into the theory. We 
have shown that since local fields contribute significantly 
to the optical and screening properties of semiconductor 
nanocrystals, the use of an atomistic quantum mechan- 
ical framework has become absolutely necessary for the 
study of nanosized structures. The present framework is 
fiexible, allowing the study of both optical and screening 
properties. It is computationally light and it gives results 
comparable to more sophisticated ab initio methods with 
good agreement with experimental data. The method 
has been applied to the study of the optical and screen- 
ing properties of both indirect gap (silicon) and direct 
gap (InAs) nanocrystals, with results that can be use- 
ful for deeper analysis. Moreover, by using the present 
framework, we demonstrate that the continuum dielec- 
tric model, almost always used to take into account local 
fields, is not adequate to study the electronic properties 
of very small nanocrystals. 
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